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(ii) That the positive electricity within the atom is nob in an electronic 
condition, but is distributed fairly uniformly through the atom. 

III. Experiments are described on the absorption of the homogeneous 
/3-rays. It is shown that the first stage in the absorption of a pencil of homo- 
geneous /3-rays consists in the scattering of the rays. The absorption of a 
completely scattered beam of homogeneous /3-rays is shown to take place 
according to an exponential law. 

In conclusion, I have much pleasure in once more recording my best 
thanks to Prof. Sir J. J. Thomson for much inspiration and advice during the 
course of these experiments. 
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(Eeceived July 8, 1910.) 

Waves of Finite Amplitude without Dissipation. 

In the investigations which follow, we are concerned with the motion of 
an elastic fluid in one dimension, say, parallel to x. It is implied not only 
that there are no component velocities perpendicular to x, but that the motion 
is the same in any perpendicular plane, so that it is a function of x and of the 
time (t) only. If u be the velocity at any point x, p the pressure, p the 
density, X an impressed force, the dynamical equation for an inviscid fluid is 

du du __ ^ 1 dp ,«k 

dt dx p dx' 

At the same time the " equation of continuity " takes the form 

dp + d(pu) _ Q> (2 , 

dt dx 

The first step, and it was a very important one, in the treatment of waves 
of finite amplitude is due to Poisson.* Under the assumption of Boyle's 
law, p = a 2 p, he proved that for waves travelling in one direction (positive) 
the circumstances of the propagation are expressed by 

u = /{#■— (a + u)t}, (3) 

* " M6moire sur la Theorie du Son/ 5 ' Journ. de l'ilcole Poly technique, 1808, vol. 7 
p. 319. 

S /t 



248 Lord Rayleigh. [July 8, 

in which / denotes an arbitrary function. When u can be neglected in 
comparison with a, this reduces to the familiar law of undisturbed propagation 
applicable to infinitesimal waves. 

Poisson does not discuss the significance of (3) further than to show that 
the boundaries of a continuous wave, limited to a finite range along x, are 
propagated with the ordinary velocity a, and that accordingly the length 
of the wave does not alter as it advances. The meaning of (3) is that in 
general u advances with a velocity equal, not to a, but to a + u, and that 
this might be expected is very easily seen (Earnshaw). From the ordinary 
theory we know that an infinitely small disturbance is propagated with a 
certain velocity a, which velocity is relative to the parts of the medium 
undisturbed by the wave. Let us consider now the case of a wave so long 
that the variations of velocity and density are insensible for a considerable 
distance along it, and at a place where the velocity (u) is finite let us 
imagine a small secondary wave to be superposed. The velocity with which 
the secondary wave is propagated through the surrounding medium is a, but 
on account of the local motion of the medium itself the whole velocity of 
advance is a-\-u, and depends upon the part of the long wave at which the 
small wave is placed. What has been said of the secondary wave applies 
also to the parts of the long wave itself, and thus we see that after a time 
t the place where a certain velocity^ is to be found is in advance of its 
original position by a distance equal, not to at, but to (a-\-u)t, or, as we may 
express it, u is propagated with velocity {a + u). 

A closer discussion of the solution represented by Poisson's integral was 
given by Stokes,* who pointed out the difficulty which ultimately arises from 
the motion becoming discontinuous. If we draw a curve to represent the 
distribution of velocity, taking x for abscissa and u for ordinate, we may find 
the corresponding curve after the lapse of time t by the following construc- 
tion : — Through any point on the original curve draw a straight line in the 
positive direction parallel to x, and of length equal to (a + u)t, or, as we are 
concerned with the shape of the curve only, equal to id. The locus of the 
ends of these lines is the velocity-curve after a time t 

But this law of derivation cannot hold good indefinitely. The crests of 
the velocity-curve gain continually on the troughs and must at last overtake 
them. After this the curve would indicate two values of u for one value of 
x, ceasing to represent anything that could actually take place. In fact we 
are not at liberty to push the application of the integral beyond the point 
at which the velocity becomes discontinuous, or the velocity-curve has a 
vertical tangent. In order to find when this happens, let us take two 
* " On a Difficulty in the Theory of Sound," ' Phil. Mag./ November, 1848. 
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neighbouring points on any part of the curve which slopes downwards in the 
positive direction, and inquire after what time this part of the curve 
becomes vertical. If the difference of abscissae be dx, the hinder point will 
overtake the forward point in the time —dxjdu. Thus the motion, as 
determined by Poisson's integral, becomes discontinuous after a time equal 
to the reciprocal, taken positively, of the greatest negative value of dujdx. 
For example, let us suppose that 

u = U cos — {x— (a + u)t], (4) 

where U is the greatest initial velocity. When t = 0, the greatest negative 
value of dujdx is — 2 7rU/\, so that discontinuity will commence at the time 
t = X/2 ttU. 

The only kind of wave travelling in the positive direction which can 
escape ultimate discontinuity is one which has no forward slope. This is 
the case of a wave forming the transition between a larger constant value 
of u when x exceeds a certain value, and a smaller constant value when x 
falls short of a certain value. As time passes, the slope everywhere becomes 
easier. We shall see presently that this wave is a wave of rarefaction, in the 
sense that during its passage the gas passes from a greater to a less density. 

It is worthy of remark that, although we may of course conceive a wave of 
finite disturbance to exist at any moment, there is in general a limit to the 
duration of its previous independent existence. By drawing lines in the 
negative instead of in the positive direction we may trace the history of the 
velocity-curve ; and we see that as we push our inquiry further and further into 
past time the forward slopes become easier and the backward slopes steeper. 
At a time equal to the greatest positive value of dxjdu, antecedent to that 
at which the curve is first contemplated, the velocity would be discontinuous. 
The exception is now a wave of condensation, involving a passage always from 
a less to a greater density. 

When discontinuity sets in, a state of things exists to which the usual 
differential equations are inapplicable; and the subsequent progress of the 
motion has not been determined. It is probable, as suggested by Stokes, that 
some sort of reflection would ensue. In regard to this matter we must be 
careful to keep purely mathematical questions distinct from physical ones. 
We shall see later how the tendency to discontinuity may be held in check 
by forces of a dissipative character. But this has nothing directly to do with 
the mathematical problem of determining what would happen to waves of 
finite amplitude in a medium, free from viscosity, whose pressure is under all 
circumstances proportional to the density.* To suppose that the problem has 

* < Theory of Sound/ 1878, § 251. 
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no solution would seem to be tantamount to admitting an inherent contra- 
diction in the assumption, usually made in hydrodynamics, of a continuous 
fluid subject to Boyle's law. It would be strange if the necessity of a 
molecular constitution for gases could be established by such an argument. 

"With Poisson's integral (3), showing how the velocity is propagated, there 
is associated another law connecting the velocity and the density in a positive 
progressive wave. In the case of a fluid obeying Boyle's law this relation is 

u— a\ogp = const. (5) 

It does not occur explicitly in Poisson's memoir, and Earnshaw considers 
that Poisson did not discover it. Certainly it is remarkable that he omitted 
to formulate the law, but at the same time it is difficult to suppose him 
ignorant of it, seeing that it follows by simple subtraction from two of his 
equations.* A formula equivalent to (5) was given explicitly, so far as I 
know, for the first time by Airy,f who attributes it to De Morgan. 

The assumption that in a progressive wave there is a definite relation 
between u and p forms the basis of Earnshaw's investigation.^ That such a 
relation is to be expected may be shown by a line of argument analogous to 
that already employed in connection with Poisson's integral. 

Whatever may be the law of pressure as a function of density, the velocity 
of propagation of small disturbances is according to the usual theory equal to 
<\/(dp/dp), and in a positive progressive wave the relation between velocity and 
condensation (s) is 

u : s = \/(dp/dp), (6) 

where s = Sp/p. If this relation be violated at any point, a wave will emerge, 
travelling in the negative direction. Let us now picture to ourselves the case 
of a positive progressive wave in which the changes of velocity and density 
are very gradual but become important by accumulation, and let us inquire 
what condition must be satisfied in order to prevent the formation of a negative 
wave. It is clear that the answer to the question whether or not a negative 
wave will be generated at any point will depend upon the state of things in 
the immediate neighbourhood of the point, and not upon the state of things 
at a distance from it, and will therefore be determined by the criterion 
applicable to small disturbances. In applying this criterion we are to 
consider the velocities and condensations, not absolutely, but relatively to 



* Equations (1), p. 364, and (b), p. 367. 

t ' Phil. Mag./ 1849, vol. 34, p. 402. The corresponding formula for long tidal waves 
of finite amplitude was also given. 

J ' Roy. Soc. Proc.,' January 6, 1859 ; ' Phil. Trans./ 1860, p. 133. 
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those prevailing in the neighbouring parts of the medium, so that the form of 
(6) proper for the present purpose is 

whence u = a / (-— ) . -^ , (8) 

J V Vfy)/ p 

which is the relation between u and /> generally necessary for a positive 
progressive wave, as laid down by Earnshaw.* 

Earnshaw worked with the so-called Lagrangian form of the equations, in 
which the motions of particular particles are followed, and he obtained 
complete solutions for a wave progressive in one direction. In the case of 
Boyle's law the relation between velocity and density is that already given (5), 
and in the case of the adiabatic law, where 

P- = U)\ (9) 

Earnshaw finds 

£)-> =1+ (r=i>, (10) 

where a is the velocity of infinitesimal disturbances under the condition 
represented by p , p 0f viz. a 2 = ypo/po- In (9) y denotes as usual the ratio of 
the two specific heats; and in (10), applicable to a, positive progressive wave, 
the constant of integration has been so chosen that u = corresponds to p = p . 
The generalised form of Poisson's integral, appropriate when p is any 
given function of p, does not appear quite explicitly in Earnshaw's 
memoir. The line of argument already used shows that it must be 

u =f[x-{u+</(dpfdp)}f]. (11) 

In the case of a gas obeying Boyle's law, 

\/(dp/dp) = const., (12) 

and (11) reduces to Poisson's form. 

In the case of the adiabatic law, we have from (9), (10), 

4) = .teY"-- «+2ri«, <u) 



V( 



<>pi \poj 

so that (Earnshaw) 

u + \/(dp/dp) = a + J (y + 1) u. (14) 

Thus (11) assumes the form 

"=/[®-{a + i(7+lM*l (15) 

and this with (10) may be considered to constitute the solution of the 

* * Theory of Sound, 5 1878, § 251. 
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problem up to the point where discontinuity sets in. We may fall back 

upon Boyle's law by putting in (15) 7=1. 

It appears that whether the relation of pressure to density be isothermal 

or adiabatic, there is a change of type as the wave advances. There can be 

no escape from such a change unless u + ^/ (dp/dp) be constant. Using (8), 

we may deduce in this case 

\/(dp/dp) = B/p, 

B being a constant, whence 

p = A-B 2 /p (16) 

expresses the only law of pressure under which waves of finite amplitude 
can be propagated without undergoing a change of type (Earnshaw). A 
simpler derivation of (16) will be given presently. 

Earnshaw further considers the genesis of disturbance in a gas originally 
at rest by the motion of a piston (supposed to be contained in a tube), but 
some of his conclusions appear to need revision. All that is required in 
these problems is virtually contained in (8), (11). If X denote the position 
of the piston at time T, the velocity of its motion is U = dX/dT, and this 
velocity is shared by the gas in contact with it. On the positive side the 
velocity of propagation of U (equal to u) is, by (11), 

TJ + x/(dp/dp); 

so that if at time t (greater than T), U is to be found at x, we must have 

x = X+iTJ + x/idp/dp)} (*-T). (17) 

Among the problems which naturally suggest themselves would be to 
determine what happens when the piston originally at rest at X = begins 
to move at time T = with a constant velocity. But if this velocity be 
positive, the discontinuity, which immediately ensues, causes the failure of 
our equations. On the other hand, if the constant velocity be negative, 
sa y — V, the initial discontinuity disappears forthwith, and the subsequent 
motion may be traced. 

Take, for example, the case of Boyle's law. We have 

X = UT, V + ^/idp/dp) = U + a; 
so that x z=TJt + a(t-T) (t > T), 

and we have to consider where the velocities and —V are to be found at 
time t. Now U = corresponds to the range of T from — 00 to 0, so that 
x ranges from at to oc . Again, U = — V corresponds to the range of T 
from to t, so that x ranges from {a— V) t to — Yt The whole range of 
x on the positive side of the piston is now accounted for, except the interval 
from x = (a— V) t to x = at. This is occupied by the transition of velocity 
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from to — V, and we infer, from what was said in the discussion of 
Poisson's integral, that this transition must take place linearly* 

Under Boyle's law the relation between velocity and density (5) is such 
that, however fast the piston may recede (u negative), a complete vacuum 
can never be formed behind it. It is otherwise under the adiabatic law (10), 
where p = corresponds to 

u = — 2a/ (y — 1) (Earnshaw), 

It may be of interest to consider further a few examples of (17). Still 
assuming Boyle's law, let us suppose that the piston is at rest (X = 0) until 
T = 0, and then moves with uniform acceleration (#), so that (T + ) 

U = #T, X = \gP = Wj2g, 

The use of these in (17) gives 

x = ^ + (U + a) t -&£*£, (18) 

2g 2y 

showing that the relation between x and U is parabolic. In (18) we see 
that dxfdU vanishes, when t = (TJ + a) fg. Thus, if g be positive, i.e. if the 
wave be one of condensation, discontinuity sets in at the front after an 
interval, reckoned from the beginning of the motion, equal to ajg. But 
if g be negative, there is no discontinuity, and (18) remains valid for an 

indefinite time. 

In general, as in the last example (g + ), the discontinuity sets in locally 
at one point of the velocity-curve, while other parts are temporarily exempt. 
It is of interest to inquire under what law the piston must advance so as to 
generate a linear velocity-curve. For then, under the adiabatic law, which 
includes Boyle's, a velocity-curve, once linear, remains linear, and if 
discontinuity enters, it must affect the whole curve simultaneously* 

It may be worth while to pause here for a moment to inquire what law 
of pressure is implied in the permanence of the linear character of the 

velocity-curve. By (7) 

dii / [dp 



9 dp V \d P J' 

so that if u + s/idp/dp) is a linear function of % du/dlogp must also be a 
linear function of u. This requires that 

pdu/dp = Cp n , 
where C and n are constants, and the most general relation between p and p 
consistent with the requirements is 

p = A + Bpr, (19) 

where A, B, 7, are constants. The relation (19) may be regarded as a kind 
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of generalised adiabatic law ; it includes the special law (16) under which 
a velocity-curve is absolutely permanent in type. 

Supposing the motion to commence at T = 0, we have 

X = flJdT = TJT- fTdTJ, 

Jo Jo 

and hence from (17), under the supposition of Boyle's law, 

x= -\TdU~-(U + a)t~-aT; (20) 

and the question before us is so to determine T as a function of U that 
(20) may be linear in IT. From (20) when t is constant 

m= -T + t-a w = 4, say, 

where t is a constant, whence 

T = *-*«> + H<r u /« 
H being the constant of integration. But, since IT = when T = 0, this 

assumes the form 

T = T' (l_e- u /*), (21) 

T' being written for t— to. Or, if we express IT in terms of T, 

V=-alog(l-^j. (22) 

In (21), (22), T' is positive, if IT is positive ; and IT becomes infinite 
when T = T'. We must therefore regard the law as limited to values of 
T less than T'. 

From (22) we find 

X = JMT = a |(T'-T)log(l-^ + T j , (23) 

which completely expresses the motion of the piston. The corresponding 
velocity-curve at time t may be verified by means of (20), (21). It is 
expressed by 

exhibiting the linear character of the slope of velocity. Evidently the slope 
becomes vertical throughout when t = T'. 

In the above example it is not necessary to suppose the law of motion of 
the piston continued up to T = T\ On the contrary, we may imagine IT to 
increase up to some prescribed finite value and then to remain constant. 
In this case the slope expressed by (24) forms the transition between IT = 
for values of x greater than at and the finite value at which the acceleration 
of the piston stops. 
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I£ the wave be one of rarefaction we must take T' negative, say — -T". In 
this case the analogue of (22) shows that U constantly increases in 
numerical value but does not become infinite in any finite time. The 
analogue of (24) is 

U = ~~,, (25) 

representing a slope which ever grows easier as time passes. 

The problem also admits of solution when the gas follows the adiabatic 
law (9). As in (15), (20), 

x = x+{a + J( 7 +l)U}(*--T).; (26) 

and (26) is to satisfy the condition of making dx/dJJ constant. In this 

dX/dU = U . dT/dTJ, 

so that ^-{(7~l)U + 24+( 7 + l)(T-0 = const. 

On integration we obtain, under the condition that U and T vanish 
together, 

/ rp \ r .1-1 T±} 

= 1, (27) 



TV 




7 + 1 
7-1 



T' being a constant which, if positive, corresponds to U = oo, thereby 
determining U as a function of T. 



For the value of X = I XJ^T, we have 
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It may be observed that, although U is infinite when T = T', X remains 
finite. Using this in (26), we get finally, on reduction, 

x = at-i(y+l)TJ(T-t), (29) 

U = -l I |^. (30) 

If ps < at, U is a linear function of x, and (T' being positive) the slope of 
the velocity curve increases until it becomes vertical when t = T'. 

If T' is negative, the wave is one of rarefaction, and (30) applies however 
great t may be. 

By putting 7=1 we fall back from (30) to (24), and less simply from 
(28) to (23).* 

* I have since found that this problem was successfully treated by Hugoniot. 



256 Lord Kayleigh. [July 8, 

Kiemann's work* is of somewhat later date than Earnshaw's, but in one 
important respect is more general It may be convenient briefly to recall 
the principal result. 

Taking p a given function of p, say <£ (p), and putting X = 0, we have 
from (1) and (2) 

<fo + u—z= -&(p) dl ° SP , dl ° gp \-u dlogp ^ -^f. 

cct (too dx at clx ctx 

If the second of these equations be multiplied by ±\/{<j>' (p)} and be 
added to the first, we find 

*--WfW>|. §=-{.-•♦'<,»£. (3D 

where 2r —f(p)-\-u, 2s=/(p)— u, (32) 

and f(p) — \\/4> f (p) • d log p* (33) 

From these follow 

^r" = — [&— {^+\/0 / (p)} ^] 

so that r remains constant when a? and tf change in such a manner that 
clx = {% + v /^> / ( i o) }^, and $ remains constant when x and t change so that 
clx = {%— -s/ft {p)}dt. In the case of a positive progressive wave s = 0, 
whence /(p) = ££ and also r = u. The velocity with which u travels in such 
a wave is accordingly u + ^S (dp/dp), of which fact (11) is merely another 
form of statement. Kiemann's equations are more general than anything 
previously given, as not limited to a single progressive wave. 

Since Eiemann's equations do not seem to have been applied in any 
example of continuous motion, I have thought it worth while to inquire 
whether they can be satisfied when r and s are both linear functions of x, 
Boyle's law being assumed, so that in (32), (33) 

</<!>' (p) = a, f(p) = a log p. 

If we suppose 

r = Ax+ B, 5=0 x+T>, (34) 

we obtain, on substitution in (31), equations for the determination of A, B, 
C, D, as functions of the time. In the first instance we find 

A-G = 1ft, (35) 

in which to t a constant may be added, and further 

A = H+1 , 0=5zLi, (36) 

* ' Gottingen Abhandhmgen/ 1860, vol. 8. 
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when H is an arbitrary constant. Also 

B-D == -aH + L/£, (37) 

L being another arbitrary constant, and thence 

B = |a(H 2 -l)log*+^J|±l> + |M, (38) 

D = -|«(E^-l)log^+ L ^7 1 ) + iy, (39) 

with M-N = - 2aH. (40) 

If we allow the origin of x to be arbitrary, as well as that of t, we may 

write 

2r = (H + l)^ + (H 2 -l)alog£ + M, (41) 

2s = (H-l)^/* + (H 2 -l)cdog* + N, (42) 

^ = r— s — x/i~~aK, (43) 

alogp = r + s = H^/^ + (H 2 -l)alog^ + J(M + N"). (44) 

If H = +1, the logarithmic term disappears, and either r or s is constant. 
In these cases we fall back upon single progressive waves. 

If H = 0, u = xjt in (43), and (44) gives 1/p proportional to t. The 
density is thus uniform with respect to x 9 but the volume of a given mass 
grows proportionally with t. The uniform expansion occurs in such a 
manner that the gas remains unmoved at the origin of co-ordinates. Since 
in this case du/dt + udu/dx = 0, we see that every part of the gas moves 
with unaccelerated velocity. 

Waves of Permanent Regime. 

When waves are propagated in one dimension without change of type, the 
circumstances are dynamically the same as in steady motion, as appears at 
once by impressing on the system a velocity equal and opposite to that of 
wave-propagation. The problem may conveniently be considered under this 
form. 

From the general equation of continuity (2), by making dp/dt equal 
to zero, or independently, we have 

p%i = p u = m, (45) 

where m is a constant which Eankine called the mass-velocity. The 
dynamical equation (1) reduces to 

u ^+lf = X. (46) 

ax a dx 
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If X = 0, (46) may be written 

? d JL=^ u 2_,. u 2 t (47) 

where uo is the velocity corresponding to p$. 
Eliminating w, we get 

' P ^ = !V(l-^?), (48) 

determining the law of pressure under which alone it is possible for a 
stationary wave to maintain itself in fluid moving (outside the wave) with 
velocity u . From (48) 



dm 



2 



-f = u i ? p A > , (49) 

dp p 2 

or p + m 2 lp = po 4- m 2 /p , (50) 

the law found by Earnshaw. 

Since, under the adiabatic law, the relation between density and pressure 
differs from (50), we conclude that a self-maintaining stationary aerial wave 
is an impossibility, unless it be in virtue of impressed forces, or of viscosity, 
or other dissipative agencies not now regarded. 

When the changes of density concerned are small, (50) may be satisfied 
approximately; and we see from (49) that the velocity of the stream 
(outside the wave) necessary to keep the wave stationary is given by 

%> = \/{dpjdp), 

which is the same as the velocity of the wave reckoned relatively to the fluid 
at a distance. 

This way of regarding the subject shows, perhaps more clearly than any 
other, the nature of the relation between velocity and density. In a 
stationary wave-form a loss of velocity accompanies an augmented density, 
according to the principle of energy, and therefore the fluid composing the 
condensed parts of a wave moves forward more slowly than the undisturbed 
portions. Relatively to the fluid at a distance, the motion of the condensed 
parts is in the same direction as that in which the waves travel. 

By means of (46), we can find what impressed force is required in order to 
ensure a stationary wave-form when (50) is not satisfied. For example, 
if p = a 2 p, we find from (45), (46), 

X = u $i + a>^R£ = (u 2_ a 2 ) dlo g n ; (51) 

showing that an impressed force is necessary at every place where it is 
variable and unequal to a. In (51) X is the accelerating force so called. 
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The actual force operative upon the element of mass p dx is Xp dx. Thus, 
on integration, 

iXpdx = m (1 )du = m(u 2 — u{\ \ 1 > , (52) 

J J\ ^7 L ^1^2 J 

if the range of integration extend from the place where the velocity is %i\ to 
the place where it becomes equal to u 2 . The integral applied force vanishes 
if the terminal velocities are such that their geometric mean is a. We may 
apply this to the case of a velocity-curve giving a simple gradual transition 
from one constant velocity U\ to another constant velocity u 2 . Under the 
above condition the integral force vanishes, but finite forces are required 
at all points of the slope, except the particular point where u = a. 

It is of some importance to notice that although, under the condition 
mu 2 = a 2 , the applied forces contribute on the whole no momentum, yet they 
do contribute energy, positive or negative. To find the work done in unit of 
time by the forces we have 

2/m . JX/) .it.dx = u 2 2 —ui 2 —2uiu 2 log (u 2 /u{). (53) 

The better to interpret this let us suppose that u^ and u 2 are positive, and 
in the first instance that u 2 >ui, so that the fluid passes from a less to a 
greater velocity, or by (45) from a greater to a less density. In the case of 
a wave of rarefaction we have therefore to consider the sign of 

j^_l_2ylogy, (54) 

when 2/>l. It is not difficult to prove that this sign is always positive. 
When y — 1 is small, the approximate value of (54) is %(y — l) 3 , and is there- 
fore positive when y>l. Again, if we remove the positive factor y from (54) 
and then differentiate, we obtain (1 — ljy) 2 , which is positive. Hence, when 
y>l, (54) is necessarily positive. The propagation of the wave of rarefaction 
without change of type requires that the impressed forces, contributing on the 
whole no momentum, should nevertheless do work upon, i.e. communicate 
energy to, the gas. 

In like manner, if the wave be one of condensation, i.e., if the gas passes 
from a less to a greater density, the operation of the impressed forces is to 
remove energy from the gas forming the wave. It follows that although 
dissipative forces, such as those arising from viscosity, may possibly consti- 
tute a machinery capable of maintaining the type of a wave of condensation, 
in no case can they maintain the type of a wave of rarefaction. 

It is desirable to extend this argument to waves propagated under the 

adiabatic law. In general, from (45), (46), 

v du , dp 

Xp = m -f + ~f\ 

so that jXpdx = m(u 2 -~ Ui)-\-p 2 — p\ = m 2 j p 2 —m 2 j p l +p 2 —p 1 . 
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As in (50), the condition that on the whole no momentum is communicated 

is 

m 2 /p2—m 2 /pi+P2--pi = 0. (55) 

Again, 

m'KiXp.u.dx = i(W-V)+ P^ = r^ ~^V^- + -V (56) 

Jjpi P hi P 2 \P 2 P l1 

The question in which we are interested is the sign of (56). If we regard 
v (the volume of unit mass), viz., 1/p, as the ordinate, andjo as the abscissa of 
a curve, the first term on the right of (56) represents the area of the curve 
bounded by two ordinates and the axis of p, while the second term is what 
the area would be if the ordinate retained throughout the mean of the 
terminal values. 

So far the argument is general. If the relation between p and v be 
adiabatic, and p2>pi, the expression (56) is negative, since v proportional to 
p-i/y makes d 2 v/dp 2 positive.* The final pressure exceeding the initial 
pressure denotes a wave of condensation, and we conclude, as before, that 
maintenance of type in such a wave requires removal of energy from the 
wave, while in the contrary case of a wave of rarefaction additional energy 
would need to be supplied. 

The problem now under discussion is closely related to one which has given 
rise to a serious difference of opinion. In his paper of 1848 already referred 
to, Stokes considered the sudden transition from one constant velocity to 
another, and concluded that the necessary conditions for a permanent regime 
could be satisfied. Eesults equivalent to his may be deduced from (45) in 
connection with the condition (u x u 2 = a 2 ) already found from (52) to express 
that there is no change of momentum on the whole. Thus, 

m = a v /(p 2 /pi\ u 2 = ct>\/(pi/p2). (57) 

Similar conclusions were put forward by Eiemann in 1860 (loc. cit). 
Commenting on these results in the " Theory of Sound" (1878), I pointed out 
that although the conditions of mass and momentum were satisfied, the 
condition of energy was violated, and that therefore the motion was not 
possible ; and in republishing this paperf Stokes admitted the criticism, 
which had indeed already been made privately by Kelvin. On the other 
hand, Burtonj and H. Weber§ maintain, at least to some extent, the original 
view. 

* Compare Lamb's * Hydrodynamics,' § 280. 
t ' Collected Works, 5 vol. 2, p. 55. 
% < Phil. Mag./ 1893, vol. 35, p. 316. 

§ ' Die Partiellen Differentialgleichungen der Mathematischen Physik,' Braunschweig, 
1901, vol. 2, p. 496. 
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Inasmuch as they ignored the question of energy, it was natural that 
Stokes and Kiemann made no distinction between the cases where energy is 
gained or lost. As I understand, Weber abandons Kiemann's solution for the 
discontinuous wave (or bore, as it is sometimes called for brevity) of 
rarefaction, but still maintains it for the case of the bore of condensation. No 
doubt there is an important distinction between the two cases; never- 
theless, I fail to understand how a loss of energy can be admitted in a motion 
which is supposed to be subject to the isothermal or adiabatic laws, in which 
no dissipative action is contemplated. In the present paper the discussion 
proceeds upon the supposition of a gradual transition between the two 
velocities or densities. It does not appear how a solution which violates 
mechanical principles, however rapid the transition, can become valid when 
the transition is supposed to become absolutely abrupt. All that I am able 
to admit is that under these circumstances dissipative forces (such as viscosity) 
that are infinitely small may be competent to produce a finite effect. 

If we suppose that under the influence of small dissipative forces the bore 
of Stokes and Eiemann can be propagated, at least approximately, we 
naturally inquire whether it can be regarded as the complete outcome of the 
simple progressive wave with a straight velocity slope which, as we have 
found, tends after a definite interval of time to assume the character of 
a bore. It would seem that the answer must be in the negative. Taking 
Boyle's law, we recognise from (5) that in the progressive wave, just before 
the formation of the bore, the relation between the velocities and densities is 

u 2 -~Ui = a log p2~~a log pi, 
while in (57) the relation is 

U2—U1 = a\/(p2/ pi) — a\/(pi/ P2). 

The two functions of p on the right, which are independent of any common 
addition to ui and u 2 , cannot be identified (unless p 2 = pi), as we have found 
already in discussing (54). This incompatibility may be regarded as 
a confirmation of Stokes' opinion that something of the nature of reflection 
must ensue. 

Permanent Regime under the influence of Dissipative Forces. 

The first investigation to be considered under this head is a very 
remarkable one by Eankine " On the Thermodynamic Theory of Waves of 
Finite Longitudinal Disturbance/'* which (except a limited part expounded 

* ' Phil. Trans., 5 1870, vol. 160, Part II, p. 27*?. 
VOL. LXXXIV. — A. T 
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by Maxwell in his " Theory of Heat ") has been much neglected.* Con- 
duction of heat is here for the first time taken into account and although 
there are one or two serious deficiencies, not to say errors, presently to be 
noticed, the memoir marks a very definite advance. 

The first step is the establishment of an equation equivalent (when the 
wave is reduced to rest) to (45), and of Earnshaw's relation (50), in which 
equations we shall usually substitute v, the volume of unit mass, for 1/p. 
Eankine remarks that " no substance yet known fulfils the condition 
expressed by (50) between finite limits of disturbance, at a constant 
temperature, nor in a state of non-conduction of heat (called the adiabatic 
state). In order, then, that permanency of type may be possible in a wave 
of longitudinal disturbance, there, must be both change of temperature and 
conduction of heat during the disturbance." However, we shall see later 
that even under Boyle's law viscosity is competent to endow a wave with 
permanency. 

The question is, how can Earnshaw's law be satisfied ? Obviously not (in 
the absence of viscosity) if the expansions are adiabatic : but if at every 
stage the right quantity of heat is added or subtracted, the gas may be made 
to follow any prescribed law. This is the idea underlying Rankine's 
investigation. For the unit mass of a perfect gas we have, as usual, pv = B0, 
denoting absolute temperature. The condition of the gas is defined by any 
two of the three quantities p, v, 0, and the third may be expressed in terms 
of them. The relation between simultaneous variations is 

cWj = dpfp + dv/v. (58) 

In order to effect the changes specified by dp and dv, it is in general necessary 
to communicate heat to the gas. Calling the necessary quantity of heat dQ, 
we may write 

M2h + CI>- (59) 

Suppose now (a) that dp = 0. Equations (58), (59), give 

:>— >-(S)i- 

where dQ/dd (p constant) expresses the specific heat of the gas under 
constant pressure. Denoting this by C, we have 

n _ fdQ\ v 

L -[di)0- 

* I must take my share of the blame. Kankine is referred to by Lamb (' Hydro- 
dynamics/ 1906, p. 466). The body of Kan kino's memoir seems to have been composed 
without acquaintance with the writings of his predecessors ; but in a supplement he 
notices the work of Poisson, Stokes, Airy, and Earnshaw. 
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Again, suppose (b) that dv = 0. We find in a similar manner that if c denote 
the specific heat under constant volume, 

_ (dQ\p 

Thus, in general, dQ/d = Gdv/v + cdp/p. (60) 

If between (58) and (60) we eliminate dp, there results 

dQ = (G-c)E^+cd6. (61) 

E 

In (61) dQ = corresponds to adiabatic expansion, when according to 
Mayer's principle the cooling effect — odd is equal to the external work done 
by the gas p dv. Hence 

C-c = E, (62) 

and therefore 7 = - = ~ , (63) 

c — E 

a relation discovered by Eankine himself in 1850. 

Eankine then applies (60) to find what heat must be communicated in 
order that the gas may follow Earnshaw's law making dp = — m 2 dv. With 
regard to (62), it appears that 

dQ = i^Tf) ipo + m*Vo-(y+l)p}. (64) 

It will be understood that under the condition now imposed of Earnshaw's 
relation, as well as of the ordinary gas law, there remains but one 
independent variable, and that the state of the gas may be expressed in 
terms of any one of the three quantities p, v, 6. 

We have next to consider how far the necessary supply of heat defined 
by (64) can be -effected by conduction. If the initial state (distinguished by 
suffix 1) and final state (with suffix 2) be of uniformity with respect to x, the 
total quantity of heat received by the gas during its passage must be zero, or 
Je£Q = 0. Hence from (64) Eankine finds 

p + m 2 v = \ (7 + 1) (pi +p 2 ). (65) 

This is a necessary condition ; but of course there is nothing so far to show 
that it is sufficient. 

In (65) p Qi v are any corresponding values of p and v within the wave, and 
we may identify them with p lt vi* 

Thus 

m\ = i(7-l)jPi + K7 + l)jP* (66) 

* We may, of course, also identify p 0i v Q with p 2 , v 2 , 

T 2 
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The velocity to lf equal to mv h is that with which the wave advances 
relatively to the fluid in state (1). And 

v£ = mV = v 1 {i(y-l)Pi + i(y + l)i>2} (67) 

gives the square of the velocity of wave-propagation relatively to fluid (1). 
The velocity of propagation of infinitely small disturbances ( p 2 nearly equal 
to pi) is given by % 2 = yp&i ; and thus a wave of finite condensation is 
propagated faster than an infinitesimal wave, and according to (67) a wave of 
finite rarefaction would be propagated slower than an infinitesimal wave. 
Moreover, there is no limit to the velocity of a wave of condensation. 

Eankine proceeds to express the absolute temperature (6) at a point where 
the pressure is p in a wave of permanent type. By Earnshaw's law (50) in 
combination with (65) 

= pv _ p w (7 + l)(Pi+ffa)-2p , g8 v 

6o PoVo po * (7 + 1) (pi +p 2 ) — 2po ? 

and for the ratio of terminal temperatures 

2 == P2 ^ (7 + l)^i + (7-l)^2 > / 6 m 

0i pi* (y + l)p2+(y—-l)pi' 

The second fraction on the right of (69) obviously represents the ratio of 
volumes v 2 /vi t or of densities pi/p2* 

In order to justify (65), it is not necessary that the terminal states be 
states of absolute uniformity. It will suffice that the temperature be there 
stationary (d0/dx = 0), which secures that no conduction of heat takes 
place there, and a state of stationary temperature usually involves a 
stationary pressure. To make the most of (65) we must apply it to the 
smallest ranges, i.e. between consecutive places where dp/dx vanishes. 

But here a question arises which Eankine does not seem to have con- 
sidered. In order to secure the necessary transfers of heat by means of 
conduction it is an indispensable condition that the heat should pass from the 
hotter to the colder body. If maintenance of type be possible in a particular 
wave as the result of conduction, a reversal of the motion will give a wave 
whose type cannot be so maintained. We have seen reason already for the 
conclusion that a dissipative agency can serve to maintain the type only 
when the gas passes from a less to a more condensed state. If this be so, 
the application which Eankine makes to a periodic wave is evidently 
prohibited. 

According to the second law of thermodynamics, the criterion whether the 
transformation is possible as the result of dissipative action is the sign of 



(71) 
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ldQ/0. If this be negative, the transformation is not possible. From (64) 
with use of (65) 

^ Q = 2m^-l)^ 1+i?2 - 2 ^ (70) 

In (68) we may give p , any corresponding values found in the wave. 
Thus$>o lies between^ sna.dp 2i and (7>1) 

(7 + l)(p 1 +^ 2 )~2p is positive. 
Accordingly*] c2Q/0 takes the same sign as 

* y « dp(pi+p 2 --p) 

^p{(y+l)( Pl +p 2 )-2p}' 

The integral (71) is evaluated without difficulty. Dropping the factor 
1/(7 + 1), and writing m = p 2 /pi, we get 

log- + 7log^±l±fc^. (72) 

7— l + (7+l) m 

It is evident that (72) changes sign when we substitute 1/m for m. 
If we expand (72) in powers of (m— -1) we find 

r72 ^(7 2 -l)(^-l) 3 , 

the terms in (w — 1), (gj — 1) 2 , disappearing. Thus, when w — 1 is positive, 
(72) begins positive. Differentiating (72) with respect to gt, we get 

1 + 7(7-1) ,7(7 + 1) ^ 

m 7+1 + (7 — V)m 7 — l + (7 + l) m 

When (73) is reduced to a single fraction, the denominatoi? is positive, 

and the numerator is 

(/-l)(w-l) 2 . 

We infer that if w>l, (72) is always positive, and that if gj-<1, (72) is 
always negative. Hence if P2^>pi, i& if the wave be one of condensation, 
the communications of heat required are such as may arise from conduction ; 
but if the wave be one of rarefaction, its permanency can in no wise be 
attained as the result of conduction. A wave of condensation here means 
a wave such that during its progress the gas passes always from a less dense 
to a more dense state, and the most important case is when the limits are 
finite, so that the passage constitutes the transition from one uniform density 
to a greater uniform density. 

Eankine proceeds to examine more particularly under what conditions 
a wave can be permanent. " In order that a particular type of disturbance 
may be capable of permanence during its propagation, a relation must exist 
between the temperatures of the particles and their relative positions, such 
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that the conduction of heat between the particles may effect the transfers of 
heat required by the thermodynamic conditions of permanence of type." 

The equation of conduction is readily found. The heat conducted in 
unit time across a layer of the gas is represented by h ddjdx, where k is a 
coefficient of conductivity which may be a function of the condition of the 
gas, here dependent on one variable. The equation of conduction is (u = + ) 

* 

d {-, d6\ I) r, dQ dQ 

dx\ dx/ Dt dx dx 

whence, if we reckon Q from the initial condition of constant pressure jpi, 

h^ = mQ. (74) 

dx 

And from (70) 

Q = jj^fa+ri-tp)* = 2 i^(^)( M ). (75) 

Also from (50), (65), 

p + m?v = J (7 + 1) (j>! ±p 2 ), (76) 

whence = g? =^ {(y+l)(pi+p»)-2p} , (77) 

and f = (7+l)(fi^ =i£ . (78) 

dp 2m d R 

Using these, we find with regard to (62) 

dx=z— d6 dv- hdp (y + 1 )(Pi+P2)-4:p ^ , 79 , 

mQdp F mc (y + 1) (p—pi) (P2—p) 

by which is determined the distribution of pressure (and thence of density 
and temperature) along the line of propagation. 

On the supposition that h is constant, Eankine integrates (79) in terms of 
logarithms. Writing 

p - \ (pi +P2) = ?, i te -pi) = qi> 

he obtains *** - - (lZ^Kgl±£g>Ti 

lie UUldlllb -7- 7 TT 00 > 

## mc(7+l) 9.1 T 

and « = * ffcdH£l±^log £±i + 2 log(l-^)"} , (80) 

^ being measured from the place where y = 0. Mathematically the wave is 
infinitely long; but practically the transition of pressure is effected in 
a distance comparable with Jc/mc (y + 1), which may be small in terms of 
ordinary standards. It is to be observed that the general character of the 
result does not depend upon the constancy of h 

Eeverting to (79), we see that the denominator on the right is positive, and 
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that the numerator is also positive for that part of the wave where p is nearly- 
equal to £(pi +£>2). Thus, for this part of the wave at any rate, p and x 
increase together ; or, since u is positive, the gas passes to a condition of 
greater density — the wave must be one of condensation. This consideration, 
as we have seen, Kankine overlooked. And a further limitation presents 
itself : since there cannot be two pressures in one place, it is evident that 
dx/dp must not change sign. The numerator in (79) must be positive over 
its ivhole range from p\ to j? 2 > and this will nob be the case if p%\p\ exceed 
(7+l)/(3— 7), equal for common gases to 1*61. The conclusion is that the 
only kind of wave, involving a transition from one uniform pressure to 
another, which can be maintained with the aid of conduction is a wave of 
condensation, and then only when the ratio of pressures does not exceed a 
moderate value. 

The next contribution to the subject upon which I have to comment is 
contained in a long and ably written memoir by Hugoniot.* This author, 
though he covers to a great extent the same ground, makes no reference to 
Stokes, Earnshaw, Eiemann, or Eankine, and but a very slight one to Poisson — 
a circumstance which increases the difficulty of comparison. Since Hugoniot 
uses the Lagrangian form of equation, his investigation runs naturally on 
the same lines as Earnshaw's, whose general solution for a single progressive 
wave is reproduced. I have already alluded to the solution of special 
problems relating to the propagation of a wave of variable type. 

The most original part of Hugoniot's work has been supposed to be his 
treatment of discontinuous waves involving a sudden change of pressure, with 
respect to which he formulated a law often called after his name by French 
writers. But a little examination reveals that this law is precisely the same as 
that given 15 years earlier by Kankine, a fact which is the more surprising 
inasmuch as the two authors start from quite different points of view. 
Kankine's investigation, as we have seen, is expressly based upon conduction 
of heat in the gas, but Hugoniot supposes his gas to be non-conducting. 
A question of some delicacy is here involved, which will repay careful 
examination. It will be convenient to give a paraphrase of Hugoniot's 
argument.f 

This argument depends upon an application of the principle of energy to a 
region bounded by two fixed planes, including the place of discontinuity. The 
work done by the fluid as it emerges with volume v 2 against the pressure p 2 is 

* 'Journal de TEcole Polyteehnique/ 1887, 1889. 
t Compare Lamb's < Hydrodynamics/ 1906, § 280. 
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p 2 v 2 . On the whole, therefore, the external work done by the passage of the 
unit of mass is p 2 v 2 —piVi* The increase of kinetic energy of the fluid is 

i-(u 2 2 ~-ui 2 ) = %m 2 (v 2 2 —vi 2 ) = K^ + ^i) (P1—P2), 

in virtue of (50), which requires that 

Pi— p 2 -\- / rn 2 (v 1 ~- v 2 ) = 0. (81) 

The sum of these is 

P2V2 —P1V1 + i (pi —P2) (V2 + vi) = J (^2 — tfi) fe +i?i). (82) 

We have next to consider the internal energy of unit of mass in the initial 
and final states. For this purpose we suppose the gas to expand adiabatically 
from its actual volume v to an infinite volume. In this expansion the work 
done by the gas is 

/»0O 

l p< to = -£" 9 (83) 

Jv 7-1 

so that the difference of internal energy in the two states is 

P2V2 _ JW>i_ ^ ^ 

y 1 y 1 

The principle of energy requires that the sum of this and (82) be zero, 

whence 

= (p 2 -pi)(vi + V2) / 85 x 

is the relation between the pressures and volumes in the two states. The 
result thus found by Hugoniot is the same as Eankine's. From Kankine's 
equation (65) 

pi + m 2 V! =p 2 + m 2 v 2 = |(7+1) (p t +p 2 ) = |- (pi +p 2 ) + 1- m 2 (v x + v 2 ), 
it follows that m 2 = j^±Z 2 - fc£i ? (86) 

Vi + V 2 V\ — V 2 

which is identical with (85). 

The first remark that I will make is that, although Hugoniot assumes that 
the transition between the two states is sudden, there is nothing in his 
argument which requires this, all that is really necessary being that the 
regime is permanent. The next remark is that, however valid (85) may be, 
its fulfilment does not secure that the wave so defined is possible. As a 
matter of fact, a whole class of such waves is certainly impossible, and I would 
maintain, further, that a wave of the kind is never possible under the 
conditions, laid down by Hugoniot, of no viscosity or heat-conduction. 

A closer examination of the process by which (85) was obtained will show 
that while the first law of thermodynamics has been observed, the second 
law has been disregarded. The crux of the matter lies in the comparison 
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of the internal energies of the incoming and outgoing gas expressed in (84). 
If (p 2 , V2) and (p v Vi) lie upon the same adiabatic, the work corresponding to 
the passage from the one state to the other is given without ambiguity 
by (84). But in the present case the two states do not lie upon the same 
adiabatic, and the work required is deduced upon the assumption that 
nothing is involved in the passage at v = 00 from one adiabatic to the other. 
What is actually there required is the communication (positive or negative) 
of an infinitesimal quantity of heat. From the point of view of the first 
law the infinitesimal quantity of heat may be neglected, but not so from the 
point of view of the second law, since the transfer is supposed to take place 
at the zero of temperature. When heat and work are distinguished, 
infinitesimal heat at zero may have a finite value. The imaginary passage 
to infinity has the advantage of leading rapidly to the required conclusion, 
but it rather tends to obscure the real nature of the process. While all the 
other items of the account are mechanical work, the passage from one 
adiabatic to the other (which may take place at constant finite volume) is a 
question of heat as distinguished from work. If during a complete cycle 
work would be lost and corresponding heat gained, the operation is 
dissipative and there need be no contradiction if viscosity or heat-conduction 
enter, but the opposite contingency of a gain of work at the expense of 
heat is excluded in all cases. The conclusion is the same as before. 
While a wave of condensation may, perhaps, maintain a permanent regime as 
the result of dissipative agencies, a permanent wave of rarefaction is excluded. 
It is remarked by Hugoniot that even when the ratio pi/^2 is infinite, 
^2/^1 does not exceed (y + l)/(y — 1), which for common gases is equal to 
about 6. A similar remark is made by Duhem,* who discusses the whole 
question with great generality. With regard to perfect gases " lorsqu'une 
quasi-onde de choc se propage au sein d'un gas parfait, le fluide le plus 
condense est tou jours en amont de l'onde et le fluide le moins condense en 
aval." But, so far as I see, neither of these authors proves that the 
propagation is possible in any case. 

It is a question of great interest to inquire what is the influence of 
viscosity and especially whether alone, or in co-operation with heat-conduction, 
it allows a wave of condensation to acquire a permanent regime. We 
proceed to consider this question on the basis of the usual equations, although 
it must be admitted that their application to conditions which are somewhat 
extreme raises points of uncertainty. 

Eeverting to our original equations, we recognise that (45) is unaffected 

* ' Zeitschrift f. Physikal. Chem.,' 1909, vol. 69, p. 169. 
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by the inclusion of viscosity, and that the change required in (46) is repre- 
sented by writing* 

Y _ 4 d I du\ 
3p dx\ dxj 

so that (46) takes the form 

du . dp 4 d I du\ A 

whence (v = 1/p) 
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p + m 2 v — - m//, — = pi + m 2 ^ = ^2 4- m 2 ^, (87) 

the terminal states (pi, vi), {p 2} v 2 )> being of uniformity, so that dvjdx there 
vanishes. From this it appears that (81), relating to the terminal states, 
holds good equally when viscosity is regarded. 

A simple example under the head of viscosity is to suppose the tempera- 
ture maintained uniform, as by a powerful radiation, so that the gas follows 
Boyle's law, making 
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(88) 

as governing the distribution of v along the line of propagation. In a wave 
of condensation v£>v>V2 9 so that the denominator on the right of (88) is 
positive. Thus when m is positive, dvjdx is negative, as should be the case. 
On integration (jn constant) 

-^ = {vi\og(vi-~v)-v 2 \og{v—v 2 )} } (89) 

the origin of x being chosen suitably. 

The transition of volumes from v± to v 2 occupies, mathematically speaking, 
the whole range from x = — 00 to # = + 00 , but practically it may be 
very sudden. Since in (88) dxjdv never changes sign, the condition of 
permanency for a condensational wave can always be satisfied, whatever may 
be the value of the ratio Vijv 2 \ or pi/p 2 , contrasting in this respect with the 
"limitation found to be necessary on Eankine's conclusion relative to heat- 
conduction. 

* Lamb's ' Hydrodynamics,' §§ 314, 316. 

t But the limitation pointed out by Hugoniot still obtains ; otherwise, one of the 
pressures would be negative. 
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As regards the velocity of wave propagation into the rarer medium, we 

have for its square 

%^ = m 2 Vi 2 = aH\/v 2 . (90) 

Eeturning to the case where heat development and viscosity are both 
regarded, we see that in virtue of (81) Hugoniot's reasoning is still 
applicable without change, and it leads to the same final relation (85) 
as was found by Eankine when heat-conduction is alone considered. 

In endeavouring to apply Eankine's method to the more general case 
where viscosity is retained, we shall find it more convenient to treat v, or 
(1/p), rather than p, as independent variable. If, as before, dQ denotes 
the total quantity of heat received by unit mass of the gas, we have from 

(60), (62), 

, i\dQ dv , dp 

or, on elimination of p by means of (87), 



f *dv t A d I dv\\ /niX 



In (91) dQ consists of two parts, the first (dQ^), with which alone Eankine 
dealt, the heat received by conduction, and the second (dQ 2 ) the heat 
developed internally under viscosity. As regards the latter, the heat 
developed in volume v and time dt is %vfi {dujdxfdt* in which we are to 
replace dt by dx/u, and u by mv, so that 

Multiplying this by (7—I) and subtracting it from (91), we get 

As in Eankine's investigation, the whole heat received by conduction in 
passing from one uniform state vi to another uniform state v 2 must vanish. 
Hence, on integrating between these limits, and dividing out the factor 
(V2—V1), we have 

(y+l)m 2 (vi + v 2 ) = 2y(pi + m 2 Vi) = 2y (p 2 + m 2 v 2 ) 

— 7 OPi +P2) + 7 m2 (^i + ^2), 
or, as in (86), m 2 (v x + v 2 ) = 7 (p x +p 2 ), 

the same relation as was found by Eankine. Introducing it into (93), we get 

* Lamb's ' Hydrodynamics,' § 341. 
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i/^ + i(Y + l)^i~^-^) = 0, (95) 



A particular case arises when we suppose the conductivity to be zero, so 
that dQi/dx vanishes throughout. We have then 

dv 

dx 

differing from (88) only by the factor \ (y + 1). On the supposition that fju 
is constant, the solution is nearly the same as in (89), and, in fact, reduces 
to it when 7=1, which represents Boyle's law. This case of no conduction 
is thus satisfactorily disposed of. Whatever be the ratio of pressures, a wave 
of condensation is always possible. 

It should be remarked, however, that the supposition of constant fju does 
consist with the facts as known for actual gases when 7 differs from unity. 
For such gases viscosity, though independent of density, varies with 
temperature, so that fi will not be constant in (95). But since /j, is always 
positive, this complication merely affects the particular form of the integral 
and not the general conclusion as to the possibility of a permanent wave. 

In general, from (94), if we reckon Qi from the terminal state v h 

(7— l)Qi = l(y+l)m 2 (v 1 —v)(v—v 2 ) + $mfiv-jr. (96) 



The equation of conduction is the same (74) as before. And for 0, from (87), 

dv 

7 dx^ 

so that with regard to (62) the equation of conduction becomes 



=-2!? = ~ J n^-rf^ + ^-m^ + l-m/A 
K K L ^7 



k 



dv fry + 1 , . , n o 1 t 4 d f dv 



me L dx L 2 7 J ° dx\ dx, 

= ±( y j t .l)m(v 1 -v)(v-v 2 ) + itiv — . (97) 

By omitting the terms containing p we may of course fall back on 
Eankine's problem. 

Equation (97), in its general form, is much more complicated than when 
either viscosity or heat-conduction is alone regarded, in consequence of the 
occurrence of the differential coefficient of the second order. In general, 
both h and //, are functions of temperature, and therefore of v ; but, according 
to Maxwell's theory, which assumes a molecular repulsion inversely as the 
fifth power of the distance, Cfjufk is independent of temperature (as well as 
of density), and takes the value f. And it would seem that this inde- 
pendence of temperature and density is general, seeing that the ratio is of 
no dimensions, at least so long as the repulsive force can be represented by 
an inverse power of the distance.* We shall write h for the above ratio and 

* Compare ' Koy. Soc. Proc.,' 1900, vol. 66, p. 68 ; < Scientific Papers,' vol. 4, p. 453. 
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assume that for a given gas it is an absolute constant. Thus, fi being 
written for fi/m, (97) takes the form — 

*l(*'S) + *'£ {^^-H- - J*(» + D(*-.)(.-^; 

(98) 

in which 'y 2 may be regarded as the dependent variable. For v 2 we shall 

write £, and if 

U = n* dgfdw, (99) 

our equation, since it contains x only through dx, may be reduced to one of 
the first order in U and £, i.e., 

where /(f) = 3 ^ + 1} ^ + ^ 2 -f -A, (101) 

and F(f) = ^(y+l)^^-^)^-^). (102) 

If U can be found as a function of f from (100), a? follows by simple 
integration of (99). 

In considering equation (100) we may conveniently regard f as the linear 
co-ordinate of a material particle of unit mass moving in a straight line 
with velocity XL The first term, UdTJ/dtj, then represents the acceleration 
of the particle ; the second, U/(f), may be regarded as a resistance, propor- 
tional to the velocity, and at the same time variable with the position (f) ; 
and the third term on the right hand represents a force, which is also a 
function of position. If t be the time in this subsidiary problem, U = dgfdt, 
and (100) may be written 

g+/(f)f = F(a (103) 

while by (99) dx = fi dt, (104) 

If fi be constant, the substitution of fit for x in (98) is obvious. 
If we take h = 0*4, j = 1*41, (101), (102), become 

/(|) = 0-641 ^ft + A/fa - 1-900, (105) 

F(0 = 0-723 Wh-\/§WS-y/h\ (106) 

It will be observed that, over the range from fi to f 2 , F (£) is positive, 
but that the sign of /(f) is doubtful. If f has the greater terminal value 
£i, /is negative; but, when it has the smaller terminal value f 2 , the 
sign depends upon the ratio £i/|> If t&is ratio <1*21, / is negative ; 
otherwise it is positive. 
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I suppose that a complete analytical solution of our equation is not to be 
expected, and it is, indeed, hardly necessary for our purpose. "What we most 
wish to know is whether a solution is possible which satisfies the prescribed 
conditions. Among these is the requirement that U in (100) vanish at both 
limits ; and even then the manner of evanescence must be such as to secure 
that x, as determined by (99), shall be infinite at these limits. As the 
problem originally presents itself, we should have the representative particle 
travelling in the negative direction from fi to f 2 , starting with no velocity 
and arriving with no velocity. It seems simpler to consider it in a modified 
form, i.e. with the motion reversed, so that it takes place in the direction of 
the force F. There is, then, no question of the particle stopping between 
the limits and returning upon its course. We may make this change, if in 
(105) we reverse the sign of /. We consider, then, the motion of the 
particle to be in the positive direction, from f 2 to £i, with zero velocity at 
both limits, the motion between f 2 and fi being aided by the force * F, 
which itself vanishes at these limits, and being also subject to a force of the 
nature of resistance, proportional to velocity. When £i/f# does not exceed 
1*21, the force is a resistance in the ordinary sense, i.e. it opposes the motion, 
and, in any case, it has this character near (and beyond) the arrival end fi. 
But when fi/f 2 exceeds 1*21, the force becomes what we may call a counter- 
resistance, and aids the motion near the initial end f 2 . As regards F, in the 
neighbourhood of each limit it becomes a force -proportional to distance 
therefrom, repulsive near £1, and attractive near £2. Thus, when £ is nearly 
equal to f 2 , 

F(f) = 0723^|=^&(f^fe); (107) 

.and when £ is nearly equal to f i, 

F(f) = 0-723 V%9^/& (£--£). (108) 

The particle, starting from f 2 , is bound to go through to £1. If it arrives at 
£1 with zero velocity, we shall have, presumably, a solution of our problem. 
It is possible, however, that on first arrival at fi, it may pass through, and 
only settle down after a number of oscillations. To this there does not 
appear to be any objection ; but if on the return from £1 it overshoots £2, it 
.can never again return to f 1, since on the left of f 2 the sign of / is negative ; 
and then our problem has no solution. On the other hand, from the nature 
of F, it is not possible for the particle passing through f x in the positive 
direction to escape returning. 

The character of the start from £ 2 can be investigated with the aid of 
.approximate equations. Thus, in (100) we may treat /(f) as constant, say 
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2a, where a may be either positive or negative, and take, as in (107), 
F (£ ) = /3 (£ — £ 2 ), where /3 is positive, so that 

U^? + 2«U-jS(f-f a ) = 0. (109) 

If in (109) we assume 

U = X(f-£ a ), (110) 

we find that the equation is satisfied provided that 

X= -*±v/(" 2 + /3), (111) 

one value (Xi) being positive and one (X2) negative. In the present case, 
where U must be positive when £>£ 2 , Xi is to be chosen. 

The differential equation (109) can be made homogeneous, and its general 
solution,* when X is real, can be put into the form 

fU-Xi(g-& £ =C (]12) 

when C is an arbitrary constant. This solution, of course, covers the cases 
where the particle starts from £ 2 with a finite velocity (U ), and it appears 
that C = U Al ~ A2 . We might conclude from this that when U = 0, then 
C = 0, but the conclusion is not safe. If, however, U and £— £ 2 are of the 
same order of magnitude, we may write U = r (£— -£ 2 ), where r is not infinite. 
Substituting this in (112), we get 

(r-XO*!. (r-X a )-**.(f-&)*i-*. = C. (113) 

When £— £ 2 vanishes, the third factor on the left is zero, and, since the first 
and second factors are not infinite, the conclusion follows that C = 0. This 
takes us back to (110), (111), the second solution (involving X 2 ) relating to 
the case where U is negative, the motion being one of approach from the 
positive side to £ 2 . 

These conclusions may be arrived at more easily from (103), of which the 
general solution in the present case is 

£-£ 2 = Ae^ + BeM, (114) 

giving U = XiAe A i' + X 2 B<M 

From these we may deduce (115) 

U-Xi(f-f 8 ) = (X a -X 1 )B^, 
U-X a (£-f a ) = (Xi->*) AeM; 

whence (112) follows by elimination of t. For our present purpose, f — £ 2 is 
to vanish when £ = — 00, so that B = ; and (110) follows with X = X x . It 
will be remarked that (110) makes x, as determined by (99), infinite when 

* See, for example, Boole's * Differential Equations,' p. 33. 
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£ = |> The circumstances of the start from £ 2 are thus definite and 
suitable. The question is as to the arrival at £1. 

In the neighbourhood of £1 the approximate equation is 

U^ + 2*'U-/3 , (£-£ 1 ) = 0, (116) 

where a! is positive and, by (108), /3' negative. If now 

U = X' (f-£0, (117) 

the values of V are \' = -*' ± vV 2 */^) ; (118) 

so that both values, if real, are negative. On the supposition of reality, 
(112) retains its form (with £1 for £ 2 ). If the velocity of arrival (TJ ) be 
finite, C = IV 1 "*' 2 , as before ; and it might be supposed that, if U vanishes 
when £--£1 = 0, C would have to vanish or become infinite. Such a con- 
clusion would be incorrect. If in (113) we suppose \\ to be numerically the 
smaller of the two values, the third factor indeed vanishes with (£—£]) as 
before ; but the conclusion that C = is evaded if ultimately r = \\. 

The situation is most easily understood from the solution in terms of t as 
in (114), (115). Since A/i, A/2 are both negative, the condition that f— fi 
and U shall vanish together when t = 00 is satisfied, whatever may be the 
values of A and B. There are now an infinite number of possible types of 
solution, instead of only one as in the former case. And it appears that the 
two simple types included under (117) are not at all upon an equal footing. 
Except in the particular case where A = 0, the solution always tends 
ultimately to the form U = \i (£— fi), and of course it may assume this form 
throughout. All these solutions satisfy the condition as to the infinitude of 
x when U = 0. 

Whether the values of A/ be real or not, the particle must ultimately settle 
down at £1, unless it escape from the region to which the approximate 
equation applies. For in (118) the real pari of \' is always negative. 

Eeturning to (100) in its general form, let us consider the variation of U 
for a given f- as dependent upon variations in /and F. We have 

d% di; 

or ^+P8U-Q = 0, (119) 

where P and Q are supposed to be known functions of £, viz., 

P = «, Q = ^p. (120) 

The solution of the linear equation (119) is 

8U = e~i Mt (iei Mt Qdg+c). (121) 
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If Q = 0, SU has the same sign as c, so that an increment of velocity 
communicated at an)'- point remains throughout of the same sign. Again, if 
Q be throughout of one sign, SU, as dependent upon it, has the same sign. 
For example, if U be positive over the range considered, Sf positive, and SF 
negative, then SU is certainly negative. The increments Bf, SF may be local, 
vanishing over any part of the range. 

The application to the present problem is obvious. If the particle passing 
any point between ^ and fi, with velocity U, arrives at £i for the first time 
without velocity, it will still arrive at £i without velocity (it must in any 
case arrive, since F is positive), if U be diminished, or if / be increased, or if 
F be diminished, or if all these changes occur together. And in the limit, 
when £i is closely approached, the ratio of U to (^ — f) is in general the 
same. 

By use of this principle we may assure ourselves as to the possibility of a 
solution in certain cases where %i/%2 does not greatly exceed unity. We 
imagine a simplified problem which admits of analytical solution and is 
derived from the actual one by alterations which everywhere (over the range 
from %2 to £1) increase F and diminish /. If this modified problem admits of 
the solution required, a fortiori will the original problem do so. 

If we consider the curve which according to (106) represents F as a 
function of £, we see that it is concave downwards, and that F will every- 
where be increased if we substitute for the curve the two terminal tangents 
at £ 2 and £1, whose equations are given in (107), (108). The abscissa of K, 
the point of intersection, is £ = ^/ (^fi). (Fig. 1.) 





MX* 
Fig. 1. Fig. 2. 

As regards /, its value is given by (105) with sign reversed, and is 
diminished when £ is diminished. The changes will therefore be in the 
required direction if we represent / over £ 2 H by its value at % 2> viz., 

f 2 = 2* = 1-259- 0-641 s, (122) 

and from H to £1 by its value at H, viz., 

/i = 2a' = 1-900-0-641 (**+*""*), (123) 

if for brevity we write s for *J (£i/£ 2 ), so that s is the ratio of terminal 
densities in the original problem. 

VOL. LXXXIV, — A. IJ 
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As regards the first portion of the course, the solution already given (110), 
(111), determines the value of U on arrival at H. We have 

IJ = Xi£ 2 (s— 1), where \i = <sj(a 2j rf3)~™ a, 

in which- a is given by (122) while by (107) ft = 0'361(s — 1). Using these, 
we get 

U = (s-l)^[^{(0'630-0'S20s) 2 ^0'361(s^l)}^0'&30 + 0'320sl (124) 

If s = 1 + a 3 where cr is small, (124) becomes 

IT =, <t% 2 x i- = 0-584&cr a . (125) 

As regards the second portion of the course, the appropriate solution is 
provided by (117), (118), where X' is restricted to be real. And in 
accordance with our suppositions a f is given by (123) while from (108) 
j3 f = — ■0 , 361(1— s" 1 ). In choosing between the values of A/ we are at 
liberty to take that which gives the largest value of U at H consistent with 
U = at fi, viz., 

Thus at H U = (s 2 -s) & K+vV 2 + /3')}> (126) 

or approximately, in terms of a (supposed small), 

U = of 2 (0*618 + 0'033(r). (127) 

From (125), (127) we may infer that when cr is small, i.e. when s does not 
much exceed unity, the particle, starting from f 2 , arrives at H with a velocity 
small enough to admit of its being stopped on arrival at fi, even under the 
simplifying conditions that have been imposed, and therefore a fortiori under 
the actual conditions of the problem. Hence when fi does not too much 
exceed f 2 , a wave of permanent regime is possible as the result of viscosity and 
heat-conduction. 

But the range of fi/f a thus proved admissible is rather severely limited. 
The postulated reality of A/ requires that ot' 2 + fi' be positive, and this again 
requires that s<T34. For values of s greater than this the motion in the 
simplified problem would become oscillatory. Calculation shows that, for 

$■= 1*34, (124) gives 

U = (*-l) fax 0-20, 

while (126) gives TJ = (s- 1) & x 0*43 ; 

so that up to this limit the particle starting from f 2 in the simplified problem, 
and therefore also in the actual problem, would arrive at f i with zero velocity. 
Up to a ratio of densities equal to 1*34, the wave of permanent regime is 
certainly possible. 

The next step in this line of procedure will be to replace the curve 
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representing F by the broken line lyPQIri formed by three of its tangents, of 
which two are the same terminal tangents as before, while the third, PQ, may 
be taken to be the horizontal tangent parallel to £i£ 2 . (Fig. 2.) By (106) 
the point L where F is a maximum is determined by 

-v/f = t(y/bWM = Wb(*+n 

and the corresponding value of F is 

0-1875 &(s- If. 

The abscissae of the points of intersection of the horizontal tangent with 
the terminal tangents are given by (107), (108). For M 

f=i&(s + l), 

and for E* f = i&(« 2 + «). 

As for / we are to take along £ 2 M the value at £ 2 ; along MN" the value 
at M ; and along Nfi the value at N", as given by (105), with sign changed. 

For the two terminal portions the solutions are of the same form as before, 
but for the middle portion a new form is required. Making F constant in (100) 
and writing 2a" for/, we find on integration 

~-4a ,/ ^ = Flog(F~2a ,, U) + 2a ,, U + C, (128) 

where C is an arbitrary constant, to be determined so as to suit the velocity 
with which the particle arrives at M. 

As before, the limiting value of s is determined by the consideration that, 
for our purpose, the arrival at fi must not be oscillatory. "We get s = 1*633 
about, and with this value of s it is not difficult to show that the velocity of 
arrival at N" is below the value prescribed by the solution for Nfi. The 
necessary conditions are thus fulfilled and we infer that the wave of uniform 
regime is possible so long as s<l*63. Although this ratio is moderate, it 
exceeds that found admissible in Kankine's problem, which leaves viscosity out 
of account. "We there found that the greatest admissible ratio of pressures 
was 1*61, which by (68) corresponds to 1*40 for the ratio of densities. Of 
course, in Eankine's problem the solution definitely fails at this point, while 
in the present problem all that we have so far proved is that the limit 
exceeds 1*63. 

From the low limiting values of s found necessary in these two cases in 
order to secure the reality of the roots of (118), it might be inferred that the 
reality would fail for the limiting motion in the neighbourhood of ffi, when & 
had a considerable value, but such is not the case. If we use the value of/ 
from (105) appropriate to the terminal point £i itself, we have 

2*' ^ 1-900-0-641 (1 + s" 1 ) = 1-259-0-641S- 1 ; 
and /3' = -0'3615 (l-s" 1 ). 



280 Lord Rayleigh. [July 8? 

Thus, even when 5 = 00 , 

*'* + & = 0-3963-0-3615 = +0*0048; 

and the roots are always real. Hence, subsidence at £1 is not ultimately 
oscillatory, but there is nothing in this argument to exclude a finite number 
of oscillations before subsidence into the region governed by the approximate 
equation. 

The method of approximation already followed might be pushed further, 
but it seems preferable to use the general method of numerical calculation 
for the solution of differential equations as formulated by Kunge.* The 
equation is that numbered (100), in which / (whose sign is to be reversed) 
and F are given by (105), (106). If we write U' = U/£ 25 ?' = f/& our 
equation takes the form 

^ = ^(s~~ */£') (vT-1)- 1-900 + ?^ij±i). (129) 

H u * vf 

The value of s being given, it is required to trace the connection between 
XT' and £', simultaneous values being denoted respectively by a and b. 
If a receive the increment h, we have to calculate the corresponding 
increment k for b. If we call the function on the right of (129) <j> (£', U'), 
Eunge gives as a first approximation to k 

hi = <f)(a + %h,b + %<])Q.h)h, (130) 

where <f>o = <£ (a, 5). The next approximation is 

Jc = h + Uh-h), (131) 

where k 2 = \{k f + k'") and ft' = fo.h, 

k" = 0(a + M + &')^ fc'" = <Ka + M + *") A - (132) 

Having determined the new simultaneous values a + h, h-\-k, we make a 

fresh departure therefrom, and so trace out the function step by step. 

In the present application the starting point is a = 1 (i.e. % = f 2 ), & = 0, 

and we have to trace the function until g' = s 2 . The initial value of $ is 

to be found from (129) by putting £' = 1. Writing IT = £ (I' — 1)> 

we get 

0oH {1-900-0-641 (l + s)}£ -£ x 0'723(s-l) = 0, (133) 

of which the positive root is to be chosen. 

The extreme admissible value of s in our present problem is 6, and the 
first and rather elaborate calculation that I have made relates to this case. 
"From (133) <f> = 3*1591, so that taking a = 1, b = 0, A = 1, we have 
7/ = 3*1591. Calculating from (130) we find ft x = 2*255, and from (132) 

k" « 1*7076, 7c'" = 2*0772, 

* See Forsyth's " Differential Equations, 5 p. 51. 
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making k 2 =• 2*6182. Hence the correction to k h viz., ^(/^—^i), is equal to 
0121, and h = 2'376. Thus, corresponding to f ' = 2, we get U' = 2*376. The 
following are the values of IT obtained successively in this way : — 



r. 


IT. 

• 


j ! 
P. 


U'. 


«'. 


TJ'. 




i 


o-ooo 


14 


7*502 


33 


1 -2942 




2 


2-376 


18 


6-807 


34 


-8532 




3 


3-903 


22 


5-674 


35 


-4160 




4 


4-986 


26 


4-24L 


35^ 


-2043 




6 


6*380 


30 


2-610 


35| 


-1012 




10 


7*526 


32 


1-736 


36 


'0004 





The correction to h is everywhere subordinate. In the last step from 
35f to 36 no correction to ki is applied. 

There is a little difficulty in tracing by this method and with full 
accuracy the final progress to zero when £' = 36. If any doubt be left, it 
may be removed by applying the former method to the course from 35 to 36, 
using the value of / appropriate to 35 and the terminal tangent as the 
representative of the curve for F. From this it appears that even if U' at 
|' = 35 were as great as 0*7468, the moving particle could not pass £' = 36. 
The conclusion is that even in this extreme case of s = 6 the solution 
exists, and that a wave of permanent regime is possible. Further, from (99) 
we see that since II and fi' are both positive, d%jdx is positive throughout, 
and the transition from the one density to the other takes place without 
alternation. 

After what has been proved little doubt could remain but that a solution 
is possible when s has any value lower than 6. I have, however, thought 
it desirable to add a rough calculation (rough on account of the relative 
magnitude of the steps) for the case of s = 3 : — 



V 

TJ' 


1 

0*00 


2 
0-90 


3 
1-19 


.iip.— ,,..- ,—..-,..— 

5 
1-26 


7 
0-82 


8 
0-40 


81 
0*19 


9 
0-01 



It is a question of some importance to consider what is the thickness of 
the transitional layer in the waves of uniform regime which have been 
proved to be possible. Mathematically speaking, the transition occupies an 
infinite space ; but if we understand the expression to refer to a transition 
approximately complete, the thickness involved is finite, and indeed 
extremely small. Eeference to (98) shows that x is of the order pf, or p/m, 
or fi/p u, where u is the velocity of the wave* For the present purpose we 
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may take u as equal to the usual velocity of sound, i.e. 3 x 10 4 cm. per 
second. For air under ordinary conditions the value of /ju/p in C.G.S. 
measure is 0*13 ; so that x is of the order ^ x 10~ 5 cm. That the 
transitional layer is in fact extremely thin is proved by such photographs as 
those " of Boys, of the aerial wave of approximate . discontinuity which 
advances in front of a modern rifle bullet ; but that according to calculation 
this thickness should be well below the microscopic limit may well occasion 
surprise. 

Resistance to Motion through Air at High Velocities. 

According to the adiabatic law the pressures and velocities in a com- 
pressible fluid free from external force, see (47), are related by 

M^l + ^^s 2 , (134) 

in which pi, pi, U\ denote the pressure, density, and velocity at one point of 
the path ; p%> p2, ih the corresponding quantities at another point. Also 
a{ 1 = ypi/pi, so that a x is the velocity of infinitesimal disturbances in the 
condition (1). In an early paper* I suggested the application of this 
formula to bodies moving through air at high velocities. Eegarding the 
obstacle as stationary and the fluid in motion with velocity u\ and pressure 
p h the pressure p% corresponding to the loss of this velocity is given by 
putting u 2 = in (133), a-\ being the ordinary velocity of sound. This is the 
pressure which should obtain at the axial point on the nose of a symmetrical 
bullet, and although this value in strictness represents the maximum 
pressure, the analogy of an incompressible fluid suggests that the mean 
pressure on a flat surface would not be greatly inferior. But in a recent 
discussion,f Mr. Mallock has shown that this formula immensely over- 
estimates the resistance actually experienced by a bullet, and (so far as I am 
aware) the discrepancy remains unexplained. 

If indeed the adiabatic law really prevailed throughout, there could be no 
escape from the conclusion formulated. A consideration of the photographs 
by Boys£ will suggest the required explanation. At a short constant 
distance in front of the bullet there is an aerial bore, or place of approxi- 
mate discontinuity, Along the axis, the fluid moving up to the bullet 
changes its density, and therefore pressure and temperature, suddenly, so 
that there is here a special opportunity for viscosity and heat-conduction 
to take effect. The pressures and velocities on the two sides of the bore are 

* < Phil. Mag,,' 1876, vol. 2, p> 430 ; < Scientific Papers,' vol. 1, p. 289. 
f <Koy. Soc, Proc./ A, 1907, vol. 79, p. 266. 

J ' Nature/ 1893, vol. 47, p. 440. The particular photograph reproduced by Mallock 
does not exhibit well the feature in question, 
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related, not according to the adiabatic law, but according to Kankine's law 
already discussed. The changes which occur may be separated into two 
stages. The first is the sudden one in which the fluid passes from the 
atmospheric condition po, p Q , with velocity u to the condition denoted by 
Pu Ph u \> After passing the bore the fluid changes gradually according to 
the adiabatic law already stated until at the nose of the bullet the condition 
is represented by p& P2> with u 2 = 0. 

We are now in a position to calculate the final pressure p$. For the first 
stage we have Kankine's formula (67), making 

or, if a? = ypo/po, & = M- ^ _ X=l , (135) 

po 7-f- 1 or 7+I 

determining the pressure just inside the bore in terms of % (the velocity of 
the bullet through quiescent air) and a, the ordinary velocity of sound. 
When uq = a, pi = p . Eor values of u less than a, the first stage does not 
exist and we may suppose p\ = po, U\ = u . 

In the second stage we use (134) with u 2 = 0. Thus 

,p\! 2 7^?i 

in which, by a formula analogous to (66), 

piui 2 = i(7+ 1 )jPo + i(7 — 1)jPi- 

Hence (a)^=<^fa^Yl + ^^V (137) 

\^o/ 47 V^o/ L 7+i^iJ 

■ When wo> a,pi/p is to be calculated from (135). When the resulting 

value is substituted in (137), P2/PQ is determined. 
When u < a, we have simply 

If ^o/a be smaZ£, (138) reduces to 



2 



Pp _ 1+ eov 

or p 2 -~ Po = 4po^o 2 , (139) 

as for an incompressible fluid. 
When u Q = &, both systems give 
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When Uofa is large, the second terms on the right of (135) and (137) may 
be neglected, and we obtain 

P3 = y+l^Hr+^\^ ; (141) 

po 2 a 2 I 47 J 

or, when we put 7 = 1*4:1, 

p 2 /p = l-30uo*/a*. (142) 

The following are some corresponding values of p 2 /po and u /a, calculated 
from (135), (137), with 7 = 1'41 :— 



u /a ...... 

JPs/Po 


1 

1-90 


2 
4-49 


3 
11-7 


4 
20-7 



From this point onwards the approximate formula (142) may suffice. The 
values found are in good agreement with Mr. Mallock's curve. 

The question as to the linear interval between the bore and the nose of the 
bullet cannot be answered from the results of the present paper. In strictly 
one-dimensional motion, the bore and the plane wall constituting the obstacle 
eould not move at the same speed. 



